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Introduction 



The solar corona e xpands non-u niformly into space as a supersonic plasma outflow 
known as the solar wind (jParkerlll958l ). The solar wind carries signatures of coronal dynam- 
ics as well as locally generated turbulent phenomena, which span a broad range of scales. 
In situ spacecraft observati ons of fluctuations in sol ar wind parameters such as velocity and 



magnet ic field (for example, iRuzmaikin et al.l (119931 ) in the ecliptic plane and lHorbury et al. 



( Il995al ) in polar fiows) typically reveal an inertial range (IR) of turbulence with a "5/3" 
inverse pow er-law scaling at high frequenc ies and a flatter "l//"-like scaling range at lower 



frequencies flMatthaeus &: Goldstein 



to evolve radially (iBavassano et al 



19861 ) . The breakpoint betw een these two ranges is seen 



19821 : iHorbury et al.l Il996al ) with the inertial range ex- 



tending to lower frequencies with increasing radial distance, suggesting a turbulent energy 
cascade in the solar wind. The solar wind also h as a background mag i ietic field and is there - 
fore a highly anisotropic plasma environment (IShebalin et al.l Il983l : lOughton et al.l Il994j ) . 
The strength of this background field relative to the amphtude of the fiuctuations deter- 
mines whether the turbulence is "strong" , i.e the amplitude of fiuctuations is comparable 
to th at of the background magnetic field (ISridhar fc GoldreichI Il994l: iGoldreich fc Sridhar 



1997 



19951) or "weak", i .e the background magnetic field is dominant (iNg fc Bhattacharjee 
Galtier et al.l 120001 ). This inertial range has been extensively studied using timeseries anal- 



(PDFs) flMarsch fc Tulll997l: IPadhve et al 



ture funtions (GSFs) (e 



20071 : IChapman fc Hnat 



ysis tech niques including power spectra (Marsch &: Tulll990a[). proba bility density functions 



200 



Horbury fc Ba 



2007 



Nicol et al 



ogh 



20081 ). 



Bruno et al. 



1997 



Hnat et al 



2004) and generalised struc- 



2005; Sorriso-Valvo et al. 



In this paper we focus on the low frequency range, where the observed a ~ 1 for 
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magnetic field fluctuations, which is ubiquitous in the solar wind and seen at all latitudes and 
radia l distances. It is dominated by signatures of coronal origin (see iMatthaeus &: Goldstein 



19861 ). unlike the turbulence seen at higher frequency, which is locally evolving. The lo- 



cation of the spectral breakpoint point between the inertial and "1//" ranges depends on 
latitude and radial distance, but it is always possible to see "1//" scaling at low frequencies. 
The power spectral density (PS D) of the "l/.f" range in the i nterplanetary magnetic field 
has been extensively studied by IMatthaeus &: Goldstein! (119861 ) and at 1 AU in the magni- 
tude of the solar wind bulk velocity v and magnetic field B by iBurlaga fc FormanI (120021 ). 



PDFs of fiuctuations in solar wind parameters at t 


lese very large scales ( \ 


Vlarsch k Tu 


1997; 


Burlaea & Forman 


2002; 


Sorriso-Valvo et al. 


2004 


; ] 


Bavassano et al. 


200. 


5) and over a 


wide 
at ions 


range of heliospheric radii. 


Burlaea & Forma 


Q ( 


200: 


I) used large scale velocity fluctu. 



at 1 AU on timescales of one hour to a year to quantify the standard deviation, kurtosis 
and skewness of PDFs over these scales. Studies of the "1//" range in the solar wind thus 
provide a unique perspective on the physics of coronal processes over the solar cycle. For 
the first time we consider components of v and B defined relative to the local magnetic 
field, and we systematically distinguish between intervals of fast and slow solar wind at solar 
maximum and minimum. Here we will focus on the anisotropy of the fiuctuations by using a 
novel decomposition technique, and will take advantage of the long timeseries available from 
the Advanced Composition Explorer {ACE) spacecraft to compare not only fast and slow 
solar wind streams but also periods of minimum and maximum solar activity. 
In the inertial range, vector quantities such as in situ velocity and magnetic field can be 
resolved for components both parallel and perpendicular with respect to the background 
magnetic field B. The duration of the ti mescale over which the background field is com- 



puted is important and both large scale B (IMatthaeus et a 



function of the scale of the fiuctuations (j Chapman fc Hnat 



19901) and average local B as a 



20071 : iHorburv et al.ll2008f ) have 



been considered in the context of inertial range turbulence. In terms of quantifying scaling, 
these approaches are generic and the focus of the present paper is to incorporate these ideas 
in statistical studies of the "1//" range, since we anticipate that coronal processes and the 
transport or propagation of coronal structures will depend on orientation with respect to 
the background magnetic field. The observed scaling would also be anticipated to depend 
quantitatively on solar cycle and to differ between fast (~ 750 km/s) or slow (^ 350km/s) 
solar wind streams. High speed fiows originate in co ronal holes ( Krieeer et al. 19731 ). whereas 
low speed fiows arise from dense coronal streamers (IGosling et al.lll98ll ). while solar rotation 
causes high and low speed fiows to interact at low latitudes. We will perform generalised 
structure function analysis (GSF) (ISornettd |2004| ) on datasets spanning these intervals in 
order to quantify the scaling properties of the magnetic and velocity field fiuctuations both 
parallel and perpendicular to the background magnetic field B. 



- 3 - 



The location of the sp ectral breakpoint be t ween the inertial and " 1/f" ranges differs in 
fast and slow streams (IHorbury et al.l l2005l : iBruno &: Carbond l2005l ) , presumably because 
at a given heliocentric distance the turbulence in the slow solar wind has had more time 
to develop than in the fast solar wind. Furthermore, the crossover between IR and "1//" 
is much clearer in fast than in slow solar wind. Here, we will see that projecting velocity 
and magnetic field parallel and perpendicular to B provides a clear indicator of where this 
crossover occurs. We compare the position of this breakpoint in fast and slow solar wind 
streams and at periods of maximum and minimum solar activity. We first see that the PSDs 
of the vector components of the velocity v and magnetic field B suggest anisotropy in the 
"1 / /" range. We then decompose v and B into parallel and perpendicular fluctuations with 
respect to the local background magnetic field B. For the simple case of quiet fast solar 
wind, we compare the PDFs of the fluctuations to see which components may or may not 
share the same underlying generating process. For completeness, we also consider the PDF 
for the density fluctuations 6p. We compare the GSFs for fast and slow solar wind at solar 
maximum and minimum. Finally, using the GSFs, we obtain values for the scaling exponents 
in the "1//" range and find that these are clearly distinct for 5f||,_L and Sb\\^±. 



2. The datasets 



The advanced composition explorer (ACE) spacecraft (jStone et al.lll998l ) orbits the La- 
grangian point sunwards of the earth (~ lAU). For the present analysis we study plasma pa- 
ramete rs (magnetic field B and velocity v) averag ed over 64 seconds from the MAG/SWEPAM 
teams (ISmith et al.lll998l : iMcComas et al.lll998l ): for the year 2007, representative of a pe- 
riod of minimum solar activity; and for the year 2000, which was a period of maximum solar 
activity. This provides datasets of ~ 4.8 x 10^ samples per year. In order to separate fast and 
slow solar wind behaviour yet still preserve a dataset with sufficient points to perform GSF 
to explore the "1//" dynamic frequency range, we divide the datasets into intervals (> 6000 
points or 4.5 days ) of fast and slow stre ams, where the cut-off between fast and slow is taken 
at 450 km/s (e.g. iHorbury et al.l 120051 ). These intervals then form one long fast solar wind 
dataset of ~ 7.4 x 10^ points, and one long slow solar wind dataset of ~ 1.4 x 10^ points 
for the year 2007 and a fast dataset of ~ 4.1 x 10^ points and a slow dataset of ~ 1.1 x 10^ 
points for the year 2000. To evaluate spectral properties, we apply Fourier techniques to the 
original continuous intervals of fast and slow solar wind. When we perform statistical anal- 
ysis using the probability density functions (PDFs) of fiuctuations in section 3, each dataset 
is treated as a single ensemble. As we preserve the time indicators for the data, the pairs of 
datapoints are always drawn from within continuous intervals of fast or slow streams. 
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We first provide an overview of tlie "1 / /" range of these data intervals by plotting the power 
spectral density F{f) of the components of v and B in the RTN coordinate system, where 
R is the sun-spacecraft axis, T is the cross product of R with the solar rotation axis, and 
is the cross product of R with T. Generally, for a signal x{t) of length A^, the power 
spectrum F{f) from the fast Fourier transform (FFT) to frequency space is given by 



for a range of frequencies / = where n = [0 : N/2] and fs is the sampling frequency. 
We take our original intervals of fast and slow solar wind and truncate (or cut) them such 
that they all have the same length of 6000 datapoints. Each interval is then split up into 
windows of 2^^ = 4096 points with a 50% overlap on the previous window. A Hamming 
window is applied to each of these sub-intervals and the FFT is computed. An average is 
then taken of these sub-interval FFTs to obtain the power spectrum for each interval. The 
power spectra for all intervals are then averaged to obtain the PSDs for fast and slow solar 
wind at both solar maximum and minimum. At lower frequencies, the magnetic field power 
spectrum F{f) ~ /° shows a spectral slope a ~ —1. Plotting F{f)/ f", a = —1 should 
therefore give a horizontal line (on average). These plots are known as compensated power 
spectra and are shown for the various solar wind conditions in Figure [TJ 




(1) 



t=i 
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Fig. 1. — Compensated power spectral density F{f)/f", a = —1 for velocity and magnetic 
field fluctuation components in the RTN coordinate system for the frequency range 10~^'^ — 
10~^Hz. Results for the fast (continous line) and slow (dashed line) are displayed separately. 
The dotted vertical lines delimit the frequency range 10~^ — 10~^Hz; this is expected to 
lie within the "1/ f" range, wi t h the breakpoint bet ween the inertial and "1//" ranges 
~ 10~^Hz ( iMarsch fc Tulll990bl : iHorbury et al.lll996al ). The three panels on the left-hand 
side are for solar maximum, while the right-hand side is solar minimum. The errors are found 
by considering one standard deviation of the datasets over which the averages are taken. 
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Figure [T] covers the expected region of transition in the spectral index of v and B 
between the IR and "1 / /" frequency ranges. However it is difficult to tell precisely whether, 
for example, the PSD behaviour between 10~^Hz and 10~^Hz really is "1//", a = 1", 
particularly in the slow solar wind. It also evident from Figure [T] that in some cases in the 
"1//"" range a varies with the solar cycle and with solar wind speed, and that for both v 
and B the a can vary from one component to another, and between v and B. This implies 
anisotropy in the fluctuations and distinct scaling between v and B. 

From a statistical point of view, let us now characterize this anisotropy by decomposing 
the velocity (or magnetic) field fiuctuations into parallel and perpendicula r components 



relative to the background magnetic field. We adopt the Taylor hypothesis (lTaylorlll938l ) 
to relate spatial and temporal scales and fiuctuations over a time lag r in the velocity (or 
magnetic field) vector components, defined as 5v{t,T) = v{t + r) — v(t). A vector average 
for the magnetic field direction 6(t, r) = B/\B\ is formed from a vector sum B[t) of all 
the observed vector B values between t — r/2 and t + 3r/2. It follows that in computing 
fluctuations over r, the background field is averaged over r' = 2r, whi ch then defines the 



minim um (Nyquist) interval necessary to capture wavelike fiuctuations (j Chapman fc Hnat 



20071 ). Using this definition of b, the inner product 

5f II = 6v h = 6vRbji + SvTbr + SvNbN (2) 

vanishes for fiuctuations which generate a velocity displacement that is purely perpendicular 
to the background magnetic field B as defined. The perpendicular fiuctuation amplitude is 
then obtained from 



Sv± = ^ 6v ■ 6v — \^6v ■ bj (3) 

We use these definitions to construct differenced timeseries 6v±(t,T), Sb±(t,T), (5f||(t,r) and 
Sb\\{t, t) over a range of r intervals within the "1//" range, that is r from a few hours up to 
a day. 



3. PDF analysis 



We first examine the PDFs of these fiuctuations and explore their possible functional 
forms. For a self-affine process, knowledge of the functional form of the PDF, and of the 
Hurst exponent H, i s sufficient in principle to build a stochastic differential e quation model 



for the process (e.g. ISornettd l2004l : IChapman et al.ll2005l : iKivaiii et al.l 120071). To comp are 



their functional form, the PDFs can be renormalised using (e.g 

P[{y- <y>)] = a-'P[a-\y~ < y >)] 



Greenhough et all 120021 ) 



(4) 
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where < ■ ■ ■ > denotes the ensemble mean and a is the standard deviation of the distribution 
From a statistical point of view, where fluctuations arise from a single physical process 
rescaling of PDFs using equation H] leads to the "collapse" of the PDFs for the different r 
onto a single function that characterizes the underlying process (e.g. Greenhough et al.ll2002 



Dudson et al.ll2005l : IPendv fc ChapmanI l2006l : IPewhurst et all 12008 



HnatetaL 



200811. Let 



us apply this technique to parallel and perpendicular velocity and magnetic field fluctuations 
in the fast solar wind at solar minimum. Figure [2] shows that the PDFs for the 6v\\ and 6vj_ 
components each collapse onto single curves that are distinct from each other. The PDF 
for 6v\\ is asymmetric about 6v\\ = 0, and we have investigated this asymmetry by sorting 
the fluctuations with respect to the sign of 6vfi into 6v'^ and The resulting GSFs and 
scaling exponents display the same fractal characteristics as 6v\\ , implying that 6v^ and 6v^ 
arise from the same physical process. 
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Fig. 2. — Parallel (upper) and perpendicular (lower) velocity fluctuations 5f ||,± in the fast 
solar wind at solar minimum for the "1//" range. The left panels show the PDFs of raw 
fluctuations sampled across intervals r between 320 and 1003 minutes. The right panels 
show the same curves normalised using equation HI 



Figure [3] shows that the PDFs for 8h\\ and each collapse onto single curves that 
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are distinct from each other. The curve for 5b\\ is distinct from that for (5f y and the PDF 
has stretched exponential tails, w hich implies that these fl uctuations may originate in multi- 
plicative or fractionating process (IFrisch fc Sornettelll997l ). The curves for 6b j_ and 6vj_ look 
remarkably similar and we will explore this later. 




Fig. 3. — Parallel (upper) and perpendicular (lower) magnetic field fluctuations in the 
fast solar wind at solar minimum for the "1//" range. The left panels show the PDFs of 
raw fluctuations sampled across intervals r between 320 and 1003 minutes. The right panels 
show the same curves normalised using equation HI 



The fun ctional fornis of t hese distributions are investigated in Figure HI A Gaussian 
distribution (IWadsworthI 119971 ) 

(5) 



1 



'2ti 



approximately fits the normalised PDFs of the 5v\\ fluctuations in the "1//" range shown in 
Figure [2] with the following fitting parameters and 95% confidence bounds: yU = ± 0.003 
(mean) and a = 1±0.002 (standard deviation). Note that since we normalised the curves to fi 
and a, an exact fit would have been fi = and a = 1 here. In contrast, the normalised PDFs 
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of the 5v± fluctuations in the "1 / /" range also shown in Figure [2] are clearly n ot Gaussian. 
Here they are fitted with three different heavy-tailed distributions: gamma (IWadsworth 
1997h defined by 

/(x|a,6) = -i^x"-^e-^'/'' (6) 

with fitting parameters from maximum likelihood estimates a = 3.083 ± 0.008 and b = 
0.580 ± 0.0 02 where the erro rs are from 95% confidence bounds; the generalised extreme 
value PDF flWadsworthlll997h defined by 



f{x\k,fi,a) 



a 



-exp 



1 + k 




1 + k- 



[X 



a 



(7) 



with fitting parameters k = 0.027 ± 0.001 (shape), fi = 0.764 ± 0.001 (location) and a = 
1.324 ±0.002 (scale). The generalised extreme value (gev) distribution combines three simple 
extreme value distributions, types I, II and III, in a single form. The value of the shape 
parameter k determines the type of the distribution. In the case k —>■ 0, the distribution is 
type I, or inverse Gumbel and equation [7| simplifies to 



f{x\k,fi, 



1 



cr 



a 



exp —exp 



a 



fi) 



a 



This distribution corresponds to a maximum extreme value distribution or the limiting dis- 
tribution of samples obtained be repeatedly selecting the maximum from an ensemble of 
even ts, which in tu rn, have a distribution with finite variance, e.g. Gaussian or exponen- 
tial flSornettd |2004| ). Types II {k > 0) and III [k < 0) of the generalised extreme value 
distribution correspond respectively to F rechet and inverse Weibull distributions. Finally a 
lognormal distribution (jWadsworthI 119971 ) defined by 



f{x\fi,a) 



xcrV 27r 



-(Inx — fi) 



(9) 



is fitted with parameters fi = 0.410 ± 0.001 and a = 0.625 ± 0.001. It can be seen from 
Figure mthat either the gamma distribution or the inverse Gumbel give good fits to the PDF 
of the Sv± fiuctuations. Physically, the gamma distribution is rela ted to the PDF of waiting 
times of events generate d by a Poisson proce ss (IWadsworthlll997l ). as noted for the analysis 
of plasma turbulence by I Graves et al.l (120021 ) . 

Turning to the magnetic field. Figure H] shows that, unlike 6v\\, for 6b\\ there is a strong 
departure from the Gaussian distribution in the tails of the PDF, which are closer to stretched 
exponential, reminiscent of turbulence. The PDFs of the 6b± fiuctuations in the "1//" 
range are fitted with the same three heavy-tailed distributions as 6v±: gamma with fitting 
parameters from maximum likelihood estimates a = 3.047 ± 0.008 and b = 0.585 ± 0.002 
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where the errors are from 95% confidence bounds; generahsed extreme value with fitting 
parameters k = 0.017 ± 0.002, = 0.770 ± 0.001 and a = 1.325 ± 0.002; and lognormal 
with fitting parameters /i = 0.406 ± 0.001 and cr = 0.633 ± 0.001. There is httle difference 
between the raw and collapsed PDFs, as 6v = foy _|_ is closer to "1//" scaling. 




Fig. 4. — Parallel and perpendicular velocity and magnetic field fiuctuations 5f ||,± and Sb\\^± 
in the fast solar wind at solar minimum for the "1//" range. For clarity, plots for only 
three representative values of r are shown for each component, whereas the fitted curves are 
computed using all the r intervals between 320 and 1003 minutes. The left panels show a 
Gaussian fit to the normalised PDF curves for both 6v\\ (upper) and 6b\\ (lower) using semilog 
y axes. The right panels show the normalised curves 6v± (upper) and 6b± (lower) fitted with 
three different distributions gamma (red); Gumbel (blue); and lognormal (grey) on semilog 
y axis. For the perpendicular components, the renormalisation with /i means that the PDFs 
are shifted so that they are centered on zero, however the gamma distribution can only have 
positive arguments. It is therefore necessary to shift the PDFs by < Sv±^rnax > /o'max- 
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Fig. 5. — Perpendicular magnetic field fluctuations Sb± ("■") and 6v± ("x") in the fast solar 
wind at solar minimum for the "1 / /" range. The gamma and Gumbel distributions are used 
to fit the normalised curves for 6b± (red) and 6vj_ (blue). Only three representative values 
of r are shown for each component, whereas the fitted curves are computed using all the r 
intervals. The left panel shows the PDFs on linear axis, whereas the right panel shows the 
normalised curves on a semilog y axis. 



As we have seen, Sb± and 6v± appear to be strongly similar in their statistics and Figure 
[5] overlays the normalised PDFs for 6b± and Sv± in the fast quiet solar wind. We see that 
they are almost identical. A possible interpretation is that both sets of fluctuations have the 
same physical process at their origin. 

For c ompleteness, we al s o exa mine the ion density fluctuations Sp in the fast quiet solar wind. 
From iMatthaeus et al.l (120071 ). one might expect these to show similar scaling behaviour to 
the Sb\\^± fluctuations, however in Figure [6] we see that this is not the case. The density 
PDFs have very sharp peaks with extended tails and are asymmetric. The rescaling collapse 
works well at the centre of the PDFs, but not towards the tails. 
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[5p-<5p>]/a [5p-<5p>]/a 

Fig. 6. — Ion density fluctuations 5p in the fast solar wind at solar minimum for the "1//" 
range. The left panel shows the PDFs of raw fluctuations sampled across intervals r between 
320 to 1003 minutes normalised using equation HI The right panel shows the normalised 
curves on a semilog y axis. 

To conclude this section, let us summarise our analysis of the PDFs of fluctuations 
in the fast quiet solar wind. Figures [2] and [3] (right-hand sides) show scaling collapse for 
j_ and 56||^_L. Figure H] (left-hand side) shows that 5v\\ and 5h\\ are distinct. This is 
manifest in both a different functional form of the rescaled PDFs and different scaling of the 
moments, which we discuss next. In particular, 5h\\ is nearly symmetric and has stretched 
exponential tails, consistent with a multiplicative process, whereas 5v\\ is more asymmetric 
and is close to Gaussian. Figure [5] shows that 5vx_ and 5h± have the same PDF functional 
form and are reasonably well fitted by the gamma and Gumbel distributions with similar 
fitting parameters, suggesting a common source for the fluctuations. 

4. GSF analysis for comparison of quiet fast and slow streams 

Scaling can be quantified by computing the generalised structure functions (GSFs) of 
the fluctuations, (| 5y-r p), where (■ ■ ■) denotes ensemble averaging, m is the order of the 
moment and 5yr = y{t + r) — y{t) is the fluctuation in a signal y{t) over a time r. Assuming 
weak stationarity and a degree of self- similarity, GSFs can be related to the scale r of the 
fluctuation by a scaling exponent C("^)? when 

Sm = (I 5y^ n ~ r^('-) (10) 



whereas the PSD measures ({2) only (e.g. iMarsch fc Tijll997l : iHorbury fc Baloghlll997l ). We 
anticipate scaling for the datasets considered here, given the indication of a "1//" range in the 
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PSDs in Figure [H however power spectr a alone cannot distinguish between fractal and multi- 
fractal behaviour (jChapman et al.ll2005l ). From equation ffTOl) . the scaling exponents ({m) are 
given quantitatively by the slopes of the GSFs. Generally for perfectly self-affine processes, 
C{m) can be described by a linear equation ({m) = Hm, where H is the Hurst exponent. 
Each successful computation of a GSF at increasingly high order yields additional information 
about the nature of the PDF of fluctuations. For practical applications of the GSF analysis to 



a broad ran g e of d atasets, see for example: MHD turbulence simulations - iMerrifield et al. 



(20051. 


2006|, 


2007); 


Chapman & Hnat 



Hnat et al 



Hnat et al 



Budaev et all (120061 ): lOewhurst et al.l (120081 ): iHnat et al. 



(120051 ): 



(|2003|) : labo - 
J2008h 



ratory plasma turbulence - 
and references therein. 

We now apply these methods to the observations. Figure [7| shows the GSFs up to m = 4 for 
6v\\, 6v±, 6b\\ and 6b± for fast and slow solar wind at solar minimum. The series is differenced 
over T = n X 640 s for n = 1 to 160, that is for a range covering ten to 1706 minutes (~ 28 
hours). The finite length of the datasets considered means that the statistics calculated 
for any given single ensemble can in principle be affected by the presence of large outliers, 
which are insufficiently n umerous to b e fully sampled. We check that this does not bias our 
results via the method of iKiyani et al.l (120071 ) . which by subtracting outliers verifies whether 
calculated exponents are robust against statistical fluctuations in the outliers. The raw and 
0.4% conditioned GSFs are shown for comparison in Figure [71 For the low-order moments 
that we consider here, we see that the difference is small, so that the finite length of our 
datasets does not significantly affect our conclusions. The raw data is used for the plots of 
the probability densities of the fluctuations in section 3. 
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Fig. 7. — Comparison of scaling properties of fluctuations in the fast solar wind (upper) and 
in the slow solar wind (lower) at solar minimum in 2007. Generalised structure functions Sm 
are plotted on log-log axes versus sampling interval r for r = 10 to 1706 minutes and m = 1 
to 4. Left panels show parallel components of fluctuations in velocity (blue) and magnetic 
field (red); right panels show corresponding perpendicular components. The raw data is 
shown by whereas "•" denote data which has been conditioned by 0.4% (the difference 
in these curves quantifies finite size effects). Linear regression fits to the "1//" range over 
r = 320 to 1002 minutes are shown. The transition from the IR to the "1//" range occurs 
at ~ 178 minutes and is shown by the dashed line. 



Figure [7] is consistent with the results shown previously, namely that v and B fiuc- 
tuations exhibit very different behaviour in the "1//" range, which corresponds to large r 
intervals. A simple self-affine noise process with PSD~ 1//", a ~ 1 would on such a plot 
have C(2) — since a = 1 + C{2). If the process is fractal then ({m) = am for all 
m. Thus we see that at r > 178 min., the GSFs for Sb\\^±_ "flatten" in the ~ 1/f range . 



consistent with previou sly reported results based on the PSD (IMatthaeus &: Goldsteinlll986 



Matthaeus et al.l 120071 ). The Sv\\^± GSFs on the contrary steepen at r > 178 minutes, show- 
ing a scaling process and exponents distinct from those of 6bii ^. These are closer to a value 
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of C(2) ~ 1, which, again for a simple noise process, is consistent with PSD~ 1//^. This is 
what we have seen in the PDF curve renormahzation of the previous section: the raw 
PDFs were close to the renormalised PDFs, since the normalization is with respect to the 
first two moments 5*1 and 5*2, which for Sb\\^± vary weakly as a function of scale r. Equation 
[To] tells one that the scaling behaviour of the process is contained in the C("^) exponents, 
given by the slopes of the GSFs. We obtain these values by linear fits to the log-log GSF 
plots. 

Whilst these results confirm the " l/.f" scaling of fluctuations in -B on long timescales, re- 
ported previously by for example iMatthaeus fc Goldstein! (Il986l ). they also highlight the 
distinct scaling of v, which we will investigate next. These GSF plots of fluctuations ori- 
ented with respect to the background field also clearly show the crossover between the IR and 
"1 / /" for fast and slow solar wind. The "1//" range is much shor t er in the slow strearn s, con- 



sistent with previous observations (e.g. iBruno &: Carbond l2005l : iHorbury et al.ll2005l ). The 



minimum value of r that we will use for the following analysis can be seen to be greater than 
the breakpoint r for both velocity and magnetic field fluctuations. For 6v^\ and Sb\\ in both 
fast and slow wind, the timescale r at which the GSFs diverge is > r = 128 minutes, the 
spectrally inferred breakpoint between IR and "1//". In contrast, the divergence between 
the GSFs of 6v± and 6b± begins at a significantly shorter timescale r ~ 100 minutes. This 
is particularly apparent when one considers the higher order moments such as m = 3, 4 
in Figure [3 It is also interesting to note that although the PDFs of 6v± and 6b± in the 
fast quiet solar wind show the same functional form (Figure [5]), their GSF scalings are very 
different. This may suggest that the fluctuations 6v± and 6b ± originate in a common coronal 
source, but their subsequent development differs in the evolving and expanding solar wind. 
Figures [8] and [9] compare the GSFs for fast and slow solar wind streams at solar maximum 
(2000) and minimum (2007); for clarity only the 0.4% conditioned results are shown. 
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Fig. 8. — Comparison of scaling properties of parallel fluctuations between fast (upper) and 
slow (lower) solar wind at solar maximum (left) in 2000 and solar minimum (right) in 2007. 
Generalised structure functions Sm are plotted on log-log axes versus sampling interval r for 
r = 10 to 1706 minutes and m = 1 to 4. Parallel components of fluctuations in velocity 
(blue) and magnetic field (red) are shown. 
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Fig. 9. — Comparison of scaling properties of perpendicular fluctuations between fast (upper) 
and slow (lower) solar wind at solar maximum (left) in 2000 and solar minimum (right) in 
2007. Generalised structure functions are plotted on log-log axes versus sampling interval 
r for r = 10 to 1706 minutes and m = 1 to 4. Perpendicular components of fluctuations in 
velocity (blue) and magnetic field (red) are shown. 

Figures [S] and [HI suggest that the scaling properties of 5v±^ and 5h±^ do not change with 
solar cycle in fast solar wind. However 5v\\ does, while the solar cycle dependence of 5h\\ is 
indeterminate. The divergences at r > 10^ in the fast solar wind at solar maximum may be 
due to finite size effects: the dataset for solar maximum is shorter than for solar minimum, 
because there are fewer long continuous time periods of fast solar wind. Figures [H] and [H] 
also show that the scaling properties of all four fluctuating quantities in the slow solar wind 
differ between solar maximum and minimum, due to different scaling exponents or a different 
breakpoint location. 

Let us summarise our conclusions from the GSF analysis. First, the breakpoints between the 
scaling properties measured by GSF analysis are different between fast and slow solar wind 
streams, and between periods of maximum and minimum solar activity. These breakpoints 
do not necessarily coincide with the breakpoint between IR and "1//" ranges inferred from 



- 18 - 



spectral analysis, however as mentioned earlier, it is difficult to extract precise quantitative 
information from the power spectra plots. The IR extends to l onger timescales in sl ow solar 



wind streams and at periods of maximum solar activity (e.g. iHorbury et al.ll2005l ). These 
trends are particularly clear in the GSFs of the perpendicular components in Figure O The 
inertial range remains relatively robust for both slow and fast solar wind streams and is 
independent of solar cycle. This is to be expected if the IR is established by a turbulent 
cascade within the evolving expanding solar wind, rather than by initial conditions in the 
corona. Intriguingly, 5v_[_ and 5h^ have the same behaviour in the "1//" range for fast solar 
wind at both solar maximum and minimum. Their scaling looks similar for the slow solar 
wind, but the breakpoint moves to longer timescales at solar maximum. All four quantities 
vary between fast and slow solar wind and solar minimum and maximum. 



5. Quantifying the scaling exponents 

Let us now quantify the observed scaling by measuring the slopes of the GSFs to obtain 
estimates of the values of the scaling exponents, C("^); the robustness of the scaling will 
also be tested. In principle, values for C("^) are obtained from the gradients of the log- log 
plots of 5*^ versus r. In practice, these are affected by the fact that both the length of the 
dataset, and the range of r over which we see scaling, are finite. As a preliminary, therefore, 
we outline a method to optimise this process to obtain a good estimate of the exponents. 



5.1. p- model and Brownian walk test timeseries 

We begin by considering a simple self-affine process where Sm ~ t^^'^\ C{m) = Hm. A 
fractal (self-affine) timeseries will in principle always give the same value of H if computed 
from any region, or range of values, of the PDF of fluctuations (differences) sampled across 
a timescale r. We seek to choose the most statistically significant subset, and we do this 
by recomputing H for different regions of the PDF; if the process is fractal, we expect to 
find the same H. To probe the full range of behaviour in the PDF, including any extended 
tails, we need to test for convergence to a single value of H for a wide dynamic range of 
the PDF, for example ~ 20cr. The largest values explored by the PDF of the data are the 
least well sampled statistically. It follows that if we successively remove these outliers, we 
should see, for a fractal timeseries, rapid convergenc e to a constan t H v alue. This is shown 



in Figure [TO] top panel for a Brownian walk, see also iKiyani et al.l (120071 ) . The scaling for a 
Brownian walk with normally distributed steps demonstrates the expected behaviour for a 
fractal process without heavy tails. On the other hand a multifractal process does not return 
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a single constant value of H as one changes th e range of values over which H is computed; 
this can be seen for the multifractal p-process (IMeneveau fc Sreenivasarull987l ) in the lower 
panel of Figure [TUl A plot of the value of the exponent (here C(2)) as we succesively remove 
outliers then can distinguish fractal and multifractal processes. For processes that are fractal, 
it also provides a more precise determination of the single exponent H that characterizes the 
timeseries. The errors are obtained by combining the least squares error in the ({m) value 
fitted across the full range, with the standard deviation of the C("^) values fitted across runs 
of data points that have varying lengths, starting with a minimum length of about half the 
total fitting range length, centered on the middle of the full fitting range. This method is 
applied for all the exponent statistics throughout this paper. 
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Fig. 10. — Finite size effects on fractal and multifractal processes. The left-hand panels 
show the scaling exponents C(^) plotted as a function of moment m = 1 to 6 for different 
percentages of removed points. The right-hand panels show C(2) plotted against the per- 
centage of removed points. The two top panels show the scaling for a Brownian walk, which 
has norm ally distributed ste ps; and the two bottom panels are the results for a multifractal 
p-model (IKiyani et al.l 120071 ). 
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5.2. Fast quiet solar wind scaling 

We now quantify the scaling exponents of 5v\\ and 5v±_ fluctuations in the fast solar wind 
at solar minimum. The corresponding GSFs are plotted in the top right panels of Figures [8] 
andO We plot the exponents C("^)) which are the gradients of the fitted power laws, from 
r = 320 to 1002 minutes in Figure [TTl In the lower panels we show how the value of C(2) 
changes as outliers are successively removed. Comparing with Figure [TOl we infer that 5v\\ 
is fractal within errors and 5v± is only very weakly multifractal (almost monofractal). For 
the exponents, we obtain C||(2) close to 1, suggestive of near Gaussian behaviour and (if the 
relation a = 1 + C(2) holds) a PSD~ 1//^. In contrast, the exponent for perpendicular 
fluctuations C±(2) is close to 0.5, implying a PSD~ 1/ p^"^. 
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Fig. 11. — Scaling properties of fluctuations 5v\\ and 5v±_ in the "1//" range, r = 320 to 1002 
minutes, in the fast solar wind at solar minimum. The upper panels show the C(^) exponents 
plotted as a function of moment m = 1 to 4 for different percentages of removed points for 
5v\\ (left) and 5v±_ (right). The bottom panel shows C(2) plotted against the percentage of 
removed points for 5v\\ (blue, upper) and 5v_[_ (red, lower). 
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Figure [12] compares the scaling exponents ({2) for 6v\\ and Sv± in fast and slow solar 
wind streams at solar minimum. The corresponding GSFs are plotted in the right-hand 
pairs of panels in Figures [S] and O Fluctuations of 6v\\ in the slow solar wind appear more 
strongly multifractal than in the fast wind. For 6v± the slow solar wind displays a much 
higher exponent value for slow (C(2) ~ 0.8) than for fast (C(2) ~ 0.5) streams, reflecting the 
intrinsic differences between the fast and slow solar wind, and the coronal plasma conditions 
and magnetic field configuration at their origin. 
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Fig. 12. — Comparisons of ({2) at solar minimum in the "1//" range as a function of the 
percentage of removed points for fast (red) and slow (blue) solar wind streams for 5v\\ (left) 
and 6v± (right). The slow solar wind scaling appears to be more strongly multifractal. 

To summarize the observed results: Analysis of the scaling exponents reveals fractal 
or weakly multi-fractal (very close to monofractal) scaling in the fluctuations of velocity 
components in the fast solar wind, with very different values for Sv\\ (C(2) ~ 0.95) and 
6v± (C(2) ~ 0.5) at solar minimum. In the slow solar wind at solar minimum, the scaling 
exponent C{2) of 6v± nearly doubles from ~ 0.5 to ~ 0.8. In contrast, the scaling of Sv\\ 
remains quantitatively similar to its value in the fast wind (Figure [T2l) . i. e. ({2) ~ 0.95 but 
has a less well-defined monofractal character. Finally, if we assume a regime in which the 
PSD /~" scaling exponent a is related to C(2) by a = 1 + C(2), then we obtain for the fast 
quiet solar wind: a ~ 1 for (as expected from Figured]), a ~ 2 for 6v\\ and a ~ 1.5 for 
Sv±. 

6. Conclusions 

We have examined the scaling of the parallel and perpendicular velocity and magnetic 
field fluctuations measured in the solar wind at ~ 1 AU by ACE, which we have decomposed 
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with respect to a locally averaged background magnetic field. Power spectra, GSFs and PDF 
collapse have been used to qualify and quantify the nature of the observed scaling in the low 
frequency "1//" range. Slow and fast solar wind streams have been compared at both solar 
maximum in 2000 and solar minimum in 2007. The slow solar wind is found to be more 
multifractal and complex than the fast solar wind. 

The magnetic field fluctuations display a flattening of the GSFs for r > 178 minutes and a 
spectral index ^ 1, consisten t with ~ 1// behaviour found previously (IMatthaeus fc Goldstein 



19861 : iMatthaeus et al.l 120071 ). In contrast, the velocity fluctuations show strong anisotropy. 



with scaling behaviour distinct from that of the B field and characterized by steepening of 
the GSFs in the "1//" range (Figure [7]) consistent with ~ 1//", a ^ 1. 
For the fast quiet solar wind, 5v\\ and 5v±_ have different scaling exponents: 5v\\ exhibits frac- 
tal scaling with C^{2) ~ 0.95 ±0.02 whereas 5vx_ is weakly multifractal with ^(2) ~ 0.49 ±0.03 
(Figure [TT]) . The PDFs for these quantities also rescale relatively well. Also in the fast quiet 
solar wind, the PDF of 6v\\ is close to Gaussian, whereas 5h\\ is nearly symmetric and has 
stretched exponential tails, consistent with a multiplicative process. The rescaled PDFs for 
5v± and 5hx_ in the fast solar wind can be fitted with the same distribution function, which 
is close to gamma or inverse Gumbel (see Figure E]). However their scaling exponents re- 
vealed by GSFs differ substantially (see Figure [9]). This is consistent with a common coronal 
source for the fluctuations but a different spatiotemporal evolution out to lAU. The func- 
tional form of the PDF then constrains the mechanism that generates the fluctuations at 
the coro na, gamma having points of contact with turbulence in confined plasmas (see for 
example iGraves et al.l l2005l : iLabit et al.l 120071 . and references therein) and Gumbel, as an 
extremal process. 

The breakpoint between the inertial range and "1//" ranges differs between fast and slow so- 
lar wind streams and between periods of maximum and minimum solar activity. The inertial 
range extends to longer timescales in slow solar wind streams and at periods of maximum so- 
lar activity. The values of the inertial range scaling exponents remain unaffected by changes 
in the solar cycle (Figure [8] and Figure [9]), consistent with locally generated turbulence. 
Our results clearly show very different behaviour between the magnetic and velocity fluctu- 
ations in the "1//" range. The fractal nature of 5v\\ points to distinct physical processes 
in the corona, and to their mapping out into the solar wind. Further work would involve 
relating the fractal scaling observed at ~ lAU with fractal stirring of magnetic footpoints in 
the corona. The different scaling observed in 5v±_ points to different dynamics perpendicular 
to the background field (field line interactions?) with a possible common coronal origin the 
the 6v±_ and 5h± fluctuations, in fast quiet solar wind. 
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